Heavily Damped Motion of One-Dimensional Bose Gases in an Optical Lattice 
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We study the dynamics of strongly correlated one-dimensional Bose gases in a combined harmonic 
and optical lattice potential subjected to sudden displacement of the confining potential. Using the 
time-evolving block decimation method, we perform a first-principles quantum many-body simula- 
tion of the experiment of Fertig et al. [Phys. Rev. Lett. 94, 120403 (2005)] across different values 
of the lattice depth ranging from the superfluid to the Mott insulator regimes. We find good quan- 
titative agreement with this experiment: the damping of the dipole oscillations is significant even 
for shallow lattices, and the motion becomes overdamped with increasing lattice depth as observed. 
We show that the transition to overdamping is attributed to the decay of superfluid flow accelerated 
by quantum fluctuations, which occurs well before the emergence of Mott insulator domains. 
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Systems of ultracold bosonic atoms provide new op- 
portunities to study many-body quantum physics in low 
dimensions, where thermal and quantum phase fluctu- 
ations play a crucial role. The unprecedented degree 
of controllability attainable in cold atonr experiments 
has led to demonstrations of striking phenomena spe- 
cific to low-dimensional systems, such as thepurely one- 
dimensional (ID) Tonks-Girardeau gases [ij, y, y] and the 
Kosterlitz-Thouless transition of 2D Bose gases [J, l5[. 

An array of ID tubes of Bose gases can be created by 
imposing a strong transverse 2D optical lattice on a Bose 
condensed gas trapped in a parabolic potential. In re- 
cent experiments of this type, transport of ID Bose gases 
has been investigated in the presence of another lattice 
potential along the ID axis by displacing suddenly the 
parabolic potential [y, |7|] , and by using a moving optical 
lattice |8| . Strong inhibition of transport with increasing 
lattice depth has been observed in these cases. While ID 
transport has also been studied in condensed nratter sys- 
tems, such as liquid "^He absorbed in nanopores [9| and 
superconducting nanowires [lO| , the flexible variability of 
optical lattice parameters allows for more detailed study 
in ultracold atomic gases. Moreover, the long relaxation 
time resulting from the diluteness of atomic gases of- 
fers the prospect of exploring dynamical properties of 
the transport over a wide range of conditions. 

Fertig et al. [7| induced dipole oscillations of ID Bose 
gases for different values of the axial lattice depth. They 
observed significant damping even for shallow lattices, 
and found a transition fronr underdamping to overdamp- 
ing with increasing lattice depth. Due to strong quan- 
tum fluctuations, mean-fleld theories fail even qualita- 
tively to describe the transport properties of ID Bose 
gases, and recent theoretical analyses of this system have 
applied approaches beyond the mean-field approxima- 
tion 111, |l^, llJ, llJ, ll5| ■ Good agreement with the exper- 
iment was obtained separately for the cases of shallow 
lattices, using a truncated Wigner approximation [12l |. 



and for deep lattices, using an extended fermionization 
method [l5[. However, so far the intermediate transi- 
tion region renrains poorly understood because the ap- 
proaches used in the previous work are no longer valid 
there. 

In this Letter, we present a direct sinrulation of this ex- 
periment from the shallow-lattice superfluid (SF) to the 
deep-lattice Mott insulator (MI) regimes, and we iden- 
tify the specific physical nrcchanism governing the tran- 
sition region. To properly include the effects of quan- 
tum ffuctuations, we use the tinre-evolving block deci- 
mation (TEBD) method, which provides precise ground 
states and real-time evolution of ID quantum lattice sys- 
tems [16|. We show that the transition to overdamping 
is a manifestation of the decay of superfiuid fiow accel- 
erated by strong quantum fluctuations, rather than that 
of the emergence of the MI domains. 



We begin with the Bose-Hubbard model [17 1 
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where b, creates a boson on the j'-th site, and fij is the 
nunrber operator. The parameters J, U^ and Q are re- 
lated to independently-measurable experimental quan- 
tities. The hopping energy, J, is expressed as J = 

AiVo/Ej^f cxp (^-C^/VoJe^) En, where A = 1.397, 

B = 1.051, and C = 2.121 Q. Vq is the axial 
lattice depth and E^^ = h? /{SmdP) is the recoil en- 
ergy, where m is the atomic mass and d is the lat- 
tice spacing. The onsitc interaction, U, is given by 
U = 2{as/d)./2nVZlE^,(VQlEYif''^ Ell, where a^ is the 
s-wave scattering length and Fl is the depth of the lat- 
tice in the transverse directions. The parabolic trap 
is centered at —Xc, and its curvature, f2, is given by 



n ~ TTiw^d /2, where wt is the frequency of the exter- 
nal harmonic potentiaL In the experiment, Xc — for 
t < 0, and it is suddenly displaced to X^. = Xq at t = 
to cause dipole motions of the Bose gas. Note that the 
Bose-Hubbard model is quantitatively valid only when 
the lattice depth is sufficiently large as Vq > 2E^ |14| . 

To compare our calculations with the experiment of 
Fertig et al. [7|, we use Os = 5.31 nm, d = 405 nm, V±_ = 
30Eji, and cut = 2tt x 60 Hz. In the experiment, the 
atoms were confined in a 2D array of decoupled tight 
ID tubes and an additional periodic potential was added 
along the direction of the tubes. We focus on the central 
tube with the number of atoms iV = 81. 

To treat the dynamics at zero temperature associated 
with Eq. ([1]), we use the TEBD method, which allows us 
to compute accurately the evolution of many-body wave 
functions of ID quantum lattice systems [16[ . We choose 
the maximum number of atoms per site rimax = 5 and 
retain up to x = 200 states in the adaptively selected 
Hilbert space (for the definition of x and a detailed con- 
vergence analysis, see Ref. [18[). Setting Xc ~ 0, we first 
calculate the ground states of Eq. ([T|) via propagation in 
imaginary time. At f = the trap center Xc is suddenly 
displaced by the distance xq to cause dipole motions, and 
we simulate the dynamics of the system via propagation 
in real time. In the experiment, the initial displacement 
was fixed to be xq — 8d and Vq/E^ was changed in the 
range of < Vq/E^ < 9 [7|. Here we also study the de- 
pendence of the damping on xq, so as to reveal the effects 
of the decay of superfluid flow. 



In Figs.[TJa)-(d), we show the local densities : 



("j), 



and their fluctuations, aj = {{n^) — n|)^/^, of the ground 
states of the system for lattice depths Vo/E-n = 2, 3, 4, 
and 8. When Vq = 2£'r, the density profile is smooth, 
with the reflected parabolic shape of the confining po- 
tential that is characteristic of SF phases. As Vb/i^R is 
increased, the density profile around Uj — 1 becomes dis- 
torted and the fluctuation aj starts to take minima at 
rij = 1 as seen in Fig. [TJc) , preluding the emergence of 
the incompressible MI domains with unit filling. When 
Vo = 8£'r (Fig. [Hd)), the unit-filling MI plateaus are 
formed and its presence is convinced by the fact that the 
value of Uj in the plateaus coincides with the value of 
the fluctuation in the homogeneous system for the unit 
filling and the same value of J/U [l9|. 

In Figs, [ije)-(h), we show the cm. position Xcm ~ 
N~^d'^, j{fij) as a function of time for xq = d and 
different values of Vq/E^. Note that this displacement 
amplitude is substantially smaller than that employed 
in the experiment of Ref. [TJ]; we return to the specific 
experimental simulation below, after identifying an im- 
portant point of physics that is seen most clearly in the 
small-displacement regime. The cm. velocity Xcm(i) for 
Vo = 3£'r is also shown by a dotted line in Fig. [Kf). 
When Vq = 2£'r, the oscillation is hardly damped; this 
means that the Bose gas behaves as a superfluid as long 



as its velocity is small. As Vq/Eb, is increased, noticeable 
damping is already seen at Vq = 3£'r. At Vq = 4i5R, the 
oscillation is drastically damped while it remains in the 
underdamped regime according to the definition of un- 
derdamping mentioned below. This means that the sup- 
pression of fluctuations at rij = 1 results in significant 
reduction of the mobility of the Bose gas even before the 
formation of the MI domains. When the lattice is very 
deep (Vq = SEn), the motion is overdamped. 

Following the method used by Ref. [7|, we model the 
cm. motion of the Bose gas in our results as a damped 
harmonic oscillator, m*x = —bx-~kx, to extract a damp- 
ing rate b/bo as a function of Vq/E-h and xo/d, where m* 



is the effective mass, k 



Urj,, 



and bo = 2mujT M- We 



regard the cm. motion as underdamped if the cm. po- 
sition oscillates across the trap center, x = ~xq. In this 
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FIG. 1: (color online) Results at four values of Vo/En'- 2 (a) 
and (e), 3 (b) and (f), 4 (c) and (g), and 8 (d) and (h). In 
the left four figures, we show the local densities rij = (rij) 
(squares) and fluctuations ctj (triangles) of the ground states. 
The dashed lines in (c) and (d) represent the fluctuations 
for the unit-flUing in the homogeneous Bose-Hubbard sys- 
tem, calculated with the use of the inflnite-size version of 
the TEBD [23|. In the right four figures, we show the time 
evolution of the cm. position a::cm(i) for a small displacement 
xo = d. In (f), icm(t) is also plotted by a dotted line. 
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FIG. 2: (color online) Damping rate 6/60 versus Vq/Eb. for 
xo = d (diamonds) and xq = 8d (squares). The triangles are 



the experimental results of Ref. [3| . 



definition, for instance, the cm. motions for Vq/E^ — 
2,3, and 4 and xq = rf shown in Figs.[lje)-(g) are under- 
damped. On the other hand, the cm. motion is regarded 
as overdamped if the first minimum of the cm. position 
does not reach the trap center as seen in Fig. [Hh). To 
obtain the damping rate, we use the minimum and max- 
imum values of Xcm and icm within the first oscillation. 
For underdamped motion, where b < bo, the damping 
rate is given by 6 = 2m*\n{Bi / B2) / At, where Bi and 
B2 are the absolute values of icm at the first and sec- 
ond extrema, respectively, and At is the time interval 
between the two extrema as indicated in Fig. [IJf). For 
overdamped motion, where b > bo, the damping rate is 
determined by solving XcmiT) = A for b, where A and 
T are the values of Xcm and t at the first maximum (see 
Fig. [T]Jh)), and Xcm{t) is the solution to the overdamped 
harmonic oscillator equation. 

In Fig. [21 we show the damping rates for xo = d and 8d 
as a function of Vo/i?R, together with the experimental 
data of Ref. [7| where Xq ~ 8d. In the experiment, the 
motion becomes overdamped with increasing the lattice 
depth, starting between Vo/E^ = 2 and 3. Our TEBD 
simulations for xo = 8d are in good agreement with the 
experimental results; the onset of the ovcrdamping lies 
between Vo/Eyi = 2.2 and 2.4, where the MI domains 
are not present. This reveals that the transition to over- 
damping observed in the experiment is not due to the 
emergence of the MI domains. Notice that although the 
damping rate b/bo exhibits a jump at the transition point, 
the actual mobility of the Bose gas is gradually reduced 
and this jump is only superficial [18 1. 

The mechanism of damping in the transition region is 
suggested in Fig. [3l which shows the maximum cm. ve- 
locities v'^^ and the damping rates for different values of 
Vq/Eyi as functions of xo/d. For a given value of the lat- 
tice depth, the damping becomes more significant as v'^^ 
increases, and the motion eventually changes from under- 
damping to overdamping. This behavior is analogous to 
the breakdown of superfluid flow observed in the experi- 
ment of Ref. [8| by using a moving optical lattice, where 
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FIG. 3: (color online) Maximum cm. velocities v^^ = \Bi\ 
(a) and damping rates b/bo (b) for Vo/En = 2 (triangles), 3 
(diamonds), and 4 (squares) as functions of xo/d. 



a transition from superfluid flow to dissi pat ive flow is in- 
duced as the lattice velocity is increased |21|. Hence, this 
strong velocity dependence of the damping rate implies 
that the transition to overdamping can be attributed to 
the breakdown of superfluid flow. 

To confirm this suggestion, we show the time evolution 
of the non-condensate fraction nnon(i) = {N — Nc)/N in 
Fig. [Ifa), where the number of the condensate bosons 
A^c is defined as the largest eigenvalue of the one-body 
density matrix (bjbi) [24]. Notice that despite the one- 
dimensionality of our system, the condensate fraction 
Nc/N of the superfiuid ground states is of order one be- 
cause our system is finite. Since as a consequence of 
the breakdown of superfluid flow, excitations kick parti- 
cles out of the condensate increasing the non-condensate 
fraction, the breakdown can be identified by growth of 
the non-condensate fraction in time as done in the ex- 
periment of Ref. [8| . In Fig. ID^a) , for small displacement 
(xq = 2d) finon hardly changes as t increases; accord- 
ingly the superfluid flow is stable. The growth of finon in 
time becomes more prominent as xo/d increases, namely, 
as the damping becomes more significant. This means 
that the cm. motion is damped by the transference of 
kinetic energy into non-condensate excitations accompa- 
nied by the breakdown of superfluid flow. In Fig. [DJb), 
one sees that the growth rate dfinon/dt takes its maxi- 
mum value immediately after the cm. velocity peaks. In 
Fig. [51 we show the maximum value of the growth rate 
dnnon/dt\jjiax for different values of the lattice depth as a 
function of Vcm^ ■ The maximum growth rate is increased 
as either v^^ or Vq/E^ is increased. When xo = 8d as 
chosen in the experiment [7|, the breakdown of super- 
fluid flow is so severe that it leads to significant damping 
already in the underdamped regime. As Vq/Eb, is in- 
creased, superfluid flow breaks down even more severely 
and the cm. motion gradually changes to overdamping. 
Therefore, the transition to overdamping observed in the 
experiment is due to the breakdown of superfluid flow. 
One can experimentally corroborate this interpretation 
by measuring the non-condensate fraction and specifying 
its growth associated with the significant damping be- 
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FIG. 4: (color online) (a) Time evolution of the non- 
condensate fraction finon(i) at Vb = 2Er for different values 
of the displacement, (b) Time evolution of dfinon/dt (solid 
line) and Ucm (dashed line) for Vb ~ 2Ek and .to ~ 8d. 



underdamping to overdamping as observed in the exper- 
iment [7| because strong quantum fluctuations accelerate 
the decay of superfluid flow, which is characterized by 
prominent growth of the non-condensate fraction in time. 
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FIG. 5: (color online) Maximum growth rates dhnon / dt\ma,x 
of the non-condensate fraction for Vq/Eb. = 2 (triangles), 3 
(diamonds), and 4 (squares) as functions of v'^^. 



cause the non-condensate fraction is one of the typical 
observables in experiments [8|. 

We now discuss the effects of quantum fluctuations on 
the breakdown of superfluid flow. As seen in Fig.[5l in the 
ID Bosc gases rffinon/c^ilmax grows gradually with increas- 
ing v'^^ in contrast to 3D Bose condensates, where a 
transition from superfluid to dissipative flow occurs very 
sharply [8|. Moreover, the superfluid flow breaks down 
signiflcantly even at much lower velocity than the criti- 
cal velocity Vcr for the dynamical instability predicted by 
mean-field theories |23l| (for instance, Vcr — 0.288i?R(i//i 
when Vo = 2£'r [24|). The smearing and reduction of 
the critical velocity stem from strong quantum fluctua- 
tions spcciflc to the ID systems [8|, which allow for the 
decay of the supercurrcnt through macroscopic quantum 
tunneling [26[. Thus, the quantum fluctuations acceler- 
ate the decay of superfluid flow, leading to the heavily 
damped motion observed in the experiment. 

In conclusion, we have studied damped dipole oscilla- 
tions of one-dimensional Bose gases in an optical lattice 
by using the time-evolving block decimation method. We 
have shown that the center of mass motion changes from 
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